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ABSTRACT 

A  roll  damping,  wind  tunnel  balance  consisting  of 
journal  and  thrust  gas  lubricated  bearings  was  de¬ 

signed  for  use  in  the  von  Karman  Gas  Dynamics  Facility's 
continuous  flow  tunnels  at  the  Arnold  Engineering  Develop¬ 
ment  Center.  The  balance  was  designed  to  be  contained  in  a 
space  envelope  of  2.5  inches  in  diameter  and  8.0  inches 
long.  It  was  required  to  support  a  radial  load  of  350 
pounds  located  ±  2  inches  from  the  journal  bearing  center 
and  a  thrust  load  of  150  pounds.  In  order  to  meet  these 
specifications,  the  bearings  had  to  be  designed  with  large 
gas  film  Reynolds  numbers  where  the  assumptions  of  laminar 
flow  and  negligible  inertia  of  the  gas  used  in  lubrica¬ 
tion  theory  may  be  invalid.  An  experimental  investigation 
of  the  bearing  resulting  from  the  initial  design  showed  it 
to  have  a  turbulent  gas  film  which  caused  an  unexpected 
detrimental  effect  on  load  capacity.  A  more  thorough 
theoretical  analysis  was  conducted  to  determine  the  effects 
of  turbulent  flow  on  gas  bearing  load  capacity.  Another 
bearing  was  designed  and  fabricated  to  operate  in  the 
laminar  flow  regime  which  improved  the  operating  character¬ 
istics  and  was  in  good  agreement  with  laminar  flow  theory. 
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CHAPTER  I 

i 

INTRODUCTION 

Inferior  performance  and  sometimes  failures  of  re¬ 
entry  vehicles  (R/V)  and  sounding  rockets  because  of  the 
coupling  of  the  roll  and  nutational  frequencies  and  other 
roll  associated  phenomena  has  generated  an  increased  inter¬ 
est  in  the  dynamic  characteristics  of  these  vehicles.  Thus 
a  roll  damping  and  Magnus  force  test  mechanism  (Figure  1) 
was  designed  and  fabricated  for  use  in  the  continuous  flow 
wind  tunnels  of  the  Von  Karman  Gas  Dynamics  Facility  (VKF) , 
Arnold  Engineering  Development  Center  (AEDC)  so  that  exper¬ 
imental  investigations  could  be  performed  to  determine  these 
roll  characteristics.  A  complete  description  of  these  test 
mechanisms  can  be  found  in  Reference  [l]1. 

A  gas  bearing  was  used  as  the  balance  in  the  roll 
damping  test  mechanism  since  it  has  very  low  frictional 
torque.  The  friction  that  does  exist  does  not  vary  with  age 
and  is  not  a  strong  function  of  the  load  which  makes  a  gas 
bearing  very  desirable  for  a  wind  tunnel  balance.  Other  gas 
bearings  have  been  used  very  successfully  at  the  VKF  as  a 


Numbers  in  brackets  refer  to  similarly  numbered 
references  in  the  bibliography. 


1 


CO 


Note:  All  components  of  the  two  mechanisms  are 
common  except  the  balances  which  are 
interchangeable. 


Figure  1.  Dynamic  Roll  Balances  used  at  the  von  Kairmah  Gas  Dynamics  Facility. 
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pitch  damping  balance  (cylindrical  bearing)  and  as  a  three- 
degree-of- freedom  balance  (spherical  bearing)  and  are  de¬ 
scribed  in  References  [1],  [2],  and  [3]. 

The  gas  bearing  to  be  used  for  the  roll  damping 
balance  was  designed  to  be  contained  in  a  space  of  2.5 
inches  in  diameter  and  8  inches  in  length.  The  balance 
was  required  to  support  both  drag  and  lift  loads,  and  since 
the  center  of  pressure  of  the  model  could  not  be  located 
exactly  at  the  journal  bearing  center,  the  balance  had  to 
sustain  pitching  moment  loads.  After  considering  all  the 
expected  configurations  to  be  tested  along  with  possible 
tunnel  conditions,  it  was  determined  that  the  balance  should 
carry  a  radial  load  of  350  pounds  located  at  +  2  inches  from 
the  bearing  center  and  an  axial  load  of  150  pounds.  The 
balance  with  a  typical  model  attached  must  also  remain  dynam¬ 
ically  stable  at  rotational  speeds  up  to  5000  revolutions 
per  minute . 

This  size  and  load  specification  required  that  the 
bearing  operate  in  a  regime  where  the  viscous  lubrication 
theory  assumptions  of  laminar  flow  and  negligible  inertia 
terms  in  the  Navier  Stokes  equations  may  be  invalid.  The 
purpose  of  this  study  is  to  cover  the  design  of  this  roll 
damping,  gas  bearing  balance  with  emphasis  on  the  gas  film 
turbulent  flow  regime  which  resulted  in  the  initial  design. 
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CHAPTER  II 

THEORETICAL  ANALYSIS 

The  theory  of  gas  film  lubrication  with  external 
pressurization  has  been  greatly  refined  since  the  first  such 
gas  bearing  was  built  in  1904  by  George  Westinghouse .  A 
history  of  gas  film  lubrication  can  be  found  in  Reference  [4] . 
Much  of  the  theoretical  advancement  for  these  bearings  has 
been  made  in  the  past  decade.  A  brief  summary  of  some  of 
these  more  recent  analyses  will  be  given,  and  the  theory 
used  in  this  analysis  will  be  presented. 

Although  there  have  been  many  different  approaches 
to  the  analysis  of  gas  bearings,  the  basic  assumptions  are 
the  same  and  are  listed  below: 

1.  The  inertia  forces  of  the  gas  film  are 
negligible . 

2.  The  gas  film  is  so  thin  that  pressure,  density, 
and  viscosity  may  be  considered  constant  across 
the  film;  and  velocity  gradients  across  the  film 
are  much  larger  than  gradients  along  the  film. 

3.  The  gas  film  flow  is  laminar. 

4.  The  gas  film  is  isothermal. 
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I.  ONE -DIMENSIONAL  ANALYSIS 

The  cylindrical  journal  bearing  shown  in  Figure  2, 
will  be  used  to  illustrate  the  different  flow  models  that 
can  be  used  in  a  theoretical  analysis .  Figure  2a  shows  a 
one-dimensional-flow  (ID)  model  that  has  been  used  exten¬ 
sively  in  early  gas  film  analysis. 

Gas  Film  Analysis 

The  analysis  of  the  flow  in  this  model,  which  is  an 
exact  model  for  a  journal  bearing  under  no  load  and  an 
infinite  number  of  feeding  orifices  (line  source  feeding) , 
is  quite  straight  forward  and  is  useful  in  illustrating  the 
basic  viscous  lubrication  theory.  Based  on  the  previously 
mentioned  assumptions,  the  Navier  Stokes  equations  for  this 
case  reduce  to 

92w  1  9p 

(1) 

9y“i  y  9z 

which  can  be  integrated  to  give 

1  9P 

w  = - y  (y-h)  . 

2y  9z 

The  mass  flow  per  unit  bearing  width  can  be  determined  to 
be 

h  h3  3p 

m'  =  f  p  wdy  - - p  —  .  (2) 

z  n  12y  9z 
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Figure  2.  Journal  Bearing  Flow  Models, 
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Introducing  the  steady  flow  continuity  equation  results  in 
the  Reynolds  equation  for  one-dimensional  flow, 


h3  9p\ 

-  p  —  =  0  (3) 

12  y  dz j 


which  describes  the  gas  film  outside  the  feeding  plane. 

Using  the  isothermal  assumption  in  Equations  (2) 
and  (3)  and  integrating  yields,  respectively,  in  dimension¬ 
less  form 


P 


2 


-1 


24  y  R-T  m'* 
z 


L 

r (—  -  £>) 
D 


(4) 


and 


P 


(-  -  £>)  + 

D  >  > 


(5) 


where  C>  is  the  maximum  of  (?,Cq) • 

The  load  normal  to  the  bearing  surface  per  unit 
length  in  the  0  direction  can  be  determined  by  integrat¬ 
ing  the  pressure,  Equation  (5),  over  the  length,  L,  which 
gives 


W" 


(1  — ““)  (3  +  Pf)"| 


3(1  +  Pf) 


rl  +  2  A)! 

L 

3 


PaL  • 

(6) 


The  only  unknown  in  the  above  equations  is  P^,  the  dimen¬ 
sionless  pressure  at  the  exit  of  the  feeding  plane,  which 
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can  be  determined  by  matching  the  mass  flow  rate  in  the  gas 
film  with  the  flow  introduced  through  the  feeding  holes. 


Flow  Restrictors 

There  are  several  types  of  inlet  restrictors,  most  of 
which  are  described  in  Reference  [5] ,  that  can  be  used  to 
regulate  the  mass  flow  through  the  feeding  plane;  but  only 
orifice  type  restrictors  with  adiabatic  flow  will  be  con¬ 
sidered  in  this  analysis.  There  are  two  types  of  restric¬ 
tors,  shown  in  Figure  3,  that  fall  in  this  category.  Figure 
3a  shows  an  "orifice  compensated"  configuration  in  which  the 
restricting  area  is  formed  by  the  orifice  cross-sectional 
area,  whereas  Figure  3b  shows  an  "inherent  compensated"  con¬ 
figuration  in  which  the  restricting  area  is  the  curtain  area 
formed  by  the  edge  of  the  orifice  and  the  film  thickness. 

The  term,  inherent  compensation,  was  first  introduced  by 
Grinnel  and  Richardson  [6] .  The  orifice  compensated  re¬ 
strictor  has  better  load  characteristics  than  the  inherent 
compensated  restrictor,  but  it  is  subject  to  an  instability 
known  as  "pneumatic  hammer"  as  found  in  References  [7],  [8], 

and  [9].  The  inherent  compensated  restrictor  does  not  ex¬ 
perience  this  instability  if  there  is  no  recess  around  the 
restrictor . 

The  orifice  mass  flow  equation  can  be  written  as 


m 


2s  .A°-Ps  G  (g£,  , 

/if-  p 

*s 


(7) 
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o 

AQ  =  7T  a  «  7T  dh 

7T  a2  «  27ra  ( d  +  h) 


b.  Inherent  compensation. 

Figure  3.  Types  of  Orifice  Restrictors. 
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where 


r 


Pf 

G  (— — ) 


2y 

Y-l 


L 


Pf 

(— > 


Y+l 

‘Pf  ~ 
(— ) 

P  = 


1/2 


when 


Pf 

(— )  >  ( - ) 

Ps  y+l 


Y 

2  Y-l 


and 


Pf 

G  (— — ) 


Y+l 

2  Y-l 

Y  ( - ) 

Y+l 


1/2 


when 


2  Y-l 


(— )  1  ( - ) 

Pc  Y+l 


The  area  term,  Aq ,  in  Equation  (7)  can  be  written  as 


A0  =' 


ira  ;  orifice  compensation 
tt  dh;  inherent  compensation. 


Equaling  Mass  Flows 


Evaluating  Equation  (4)  at  the  feeding  plane  and 
equating  the  film  mass  flow  rate  to  the  orifice  mass  flow 
rate,  Equation  (7) ,  yields 


(8) 
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p:  = 


2AV (L  -  Lx)  (1  +  e  cos  6) 


D 


-a 


G(t>  + 1  ' 


(9) 


where 


6Cn  yNv'ST  „  ( a-2 )  (3-a) 

A  =  — - - —  (a  ;  (d) 


P  c 
*a 


a 


(9a) 


h  =  c  (1  +  e  cos  0)  , 


e  =  A/c  , 


(Figure  2c,  page  6) 


a  = 


'3  for  orifice  compensation 
.2  for  inherent  compensation. 


The  stiffness  at  e  =  0  can  be  determined  by  the  relation. 


dW 

de 


=  R  I 

e=0  0 


2ir  dW"  dP , 


dP^  de 


cos  0  d0  . 


(10) 


e=0 


The  first  derivative  can  be  evaluated  from  Equation  (5)  in 
terms  of  Pf  to  give 


1  dW"* 


L 


p  L  dP* 
*a  f 


5Dl=i- 


1  + 


a  +  pf)2j 


(id 


The  second  derivative  in  Equation  (10)  can  be  determined 
from  Equation  (9)  to  be 
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(12) 


where  the  term  in  braces  is  1  for  choked  flow.  Since 

and  are  independent  of  0  Equation  (10)  can  be  integrated 

to  yield  in  dimensionless  form 


ir  a 

=  -  (D  )  (D  )  (13) 

e=0  2 (V  -  1)  1  Z 

which  is  a  function  of  Pf,  Pf  can  be  related  to  A  by  using 
Equation  (9)  (e  =  0)  and  the  stiffness  can  thus  be  deter¬ 
mined  as  a  function  of  the  restrictor  coefficient.  Tang  and 
Gross  used  a  similar  analysis  and  their  design  charts  are 
presented  in  Reference  [10].  It  can  be  seen  from  Equation 
(13)  that  the  stiffness  of  a  gas  bearing  with  inherent  com¬ 
pensation  is  two-thirds  the  stiffness  of  the  same  configura¬ 
tion  with  orifice  compensation  for  a  fixed  value  of  P^. 

Although  the  one-dimensional-flow  analysis  does  pro¬ 
vide  a  reasonably  quick  method  to  obtain  rough  design  esti¬ 
mates  for  some  bearing  configurations,  it  usually  predicts 
the  optimum  mass  flow  rate  to  be  too  large  (too  large 
orifice  diameter)  and  overestimates  load  capacity  by  as 
much  as  100  percent. 


dW 

de 
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II.  POINT  SOURCE  EFFECT 

Gas  bearing  films  are  most  often  fed  through  discrete 
orifices  (point  source  feeding)  as  shown  in  Figure  2b,  page 
6,  rather  than  the  line  source  discussed  previously.  This 
is  the  prime  reason  for  the  ID  analysis  to  predict  orifice 
diameters  too  large  for  maximum  stiffness.  The  point  source 
effect  has  been  investigated  by  Dudgeon  and  Lowe  [11]  using 
an  analog  analysis.  The  results  from  this  investigation 
were  used  to  provide  corrections  to  the  line  source  analysis 
The  analysis  used  in  Reference  [6]  also  used  a  line  to  point 
source  correction,  which  was  obtained  from  data  presented 
in  Reference  [12].  These  corrections  help  the  one-dimen¬ 
sional-flow  analysis  to  give  more  accurate  predictions  but 
it  still  usually  overestimates  the  load  capacity. 

III.  CIRCUMFERENTIAL  FLOW  EFFECT 

When  a  load  is  applied  to  the  journal  bearing,  the 
gas  in  the  film  begins  to  flow  from  the  high  pressure 
(small  clearance)  side  of  the  bearing  to  the  low  pressure 
(large  clearance)  side  as  shown  in  Figure  2c,  page  6. 

This  flow  tends  to  reduce  the  circumferential  pressure 
gradient  and  thus  load  capacity  predicted  by  the  ID  analysis 
The  error  introduced  by  ignoring  the  circumferential  flow 
becomes  larger  as  the  bearing  length  to  diameter  ratio  is 
increased. 
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Robinson  and  Sterry  [13]  and  Shires  [14]  applied 
empirical  corrections  to  account  for  the  circumferential 
flow.  Dudgeon  and  Lowe  [11]  used  the  average  pressure  ob¬ 
tained  from  the  one-dimensional-flow  analysis,  applied  this 
pressure  to  obtain  the  circumferential  flow,  and  iterated 
to  obtain  the  load  capacity.  Lund  et  al .  [15]  used  an 

e-perturbation,  two-dimensional-flow  (2D)  analysis  which 
will  be  discussed  later  to  account  for  the  circumferential 
flow.  In  one  case  this  analysis  was  used  with  an  empirical 
correction  for  the  point  source  effect  to  give  very  good 
results.  Extensive  design  data  using  this  analysis  are 
available  in  Reference  [16] .  In  another  case  a  more  exact 
analysis,  which  will  be  used  in  this  investigation,  was 
obtained  using  the  e-perturbation  technique  considering 
point  source  flow. 

IV.  e-PERTURBATION,  POINT  SOURCE  JOURNAL 
BEARING  ANALYSIS 

Before  the  e-perturbation  analysis  is  applied,  it  is 
necessary  to  obtain  the  two-dimensional-flow  Reynolds  equa¬ 
tion.  Using  the  basic  gas  film  assumptions,  the  Navier 
Stokes  equation  in  the  circumferential  direction  reduces 
to 

32u  i  3p 

3y2  Ur  38 
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and  similar  to  the  ID  analysis,  the  circumferential  velocity 
and  mass  flow  rate  per  unit  bearing  length  can  be  expressed 
as 


1  9P 

u  = - y(y  -h)  , 

2y  r90 


mQ  =  - 


hJ  9  p 

-  p  - 

12y  r  9  6 


(14) 


Using  Equations  (2)  and  (14)  along  with  the  continuity  equa¬ 
tion  results  in  the  Reynolds  equation  for  two-dimensional- 
flow, 


9 

9z 


-  hJ 

9p  ~ 

P  — 

1 

+  — y 

9 

~  hJ 

9P1 

_12  y 

9z  _ 

r^ 

96 

_12y 

96  J 

J  +  =  0  , 


where  is  the  external  mass  flow  rate  per  unit  surface 

area.  Using  the  isothermal  assumption  and  introducing  the 
dimensionless  parameters. 


P  = 


h  z 

H  =  —  ,  and  £  =  —  , 
c  r 


gives 


9 

r  3  3P' i 

9 

r  3  9P  ~| 

— 

H  - 

+  — 

L  9?  J 

96 

L  96  J 

2  4r  y  R-T 

3  2  mf  =  0 

c  Pa 
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Since  in  this  analysis  the  orifices  are  represented  by  point 
sources,-  the  feeding  term  can  be  replaced  by  the  source 
strength  and  the  preceding  equation  becomes 


3 

r  3  9p2“ 

3 

r  3  9p2l 

— 

H  - 

+  - 

H  - 

35 

l  35- 

30 

L  30  J 

N 

-  6C5  -  50)  l  C_.  6(0  -  0_.)  , 
j=l 


(15) 


where 


6(  )  is  the  Dirac  delta  function, 

£q  is  the  ^-coordinate  of  the  point  sources, 

th 

0^  is  the  angular  position  of  the  j  source, 

N  is  the  number  of  sources  on  each  side  of  the 


origin 


.  th 


Cj  is  the  strength  of  the  j  point  source. 


Since  the  bearing  is  symmetric  in  £,  it  is  sufficient  to 
consider  only  half  of  the  bearing.  An  unfolded  view  of  the 
bearing  along  with  the  coordinate  system  and  boundary  con¬ 
ditions  is  shown  in  Figure  4. 

It  is  clear  from  Figure  4a  that  H  can  be  expressed 


as 


H  =  1  +  e  cos  0  (16) 

where- e  is  defined  as  A/c .  Now  expanding  P  and  in 
power  series  of  e  and  omitting  all  terms  in  e  of  secbnd 
degree  or  higher  gives 
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P2  =  +  2e  PQ  PL  ,  (17a) 

H3  =  1  +  3e  cos  9 ,  (17b) 

C.  =  C?  +  £  C*  .  (17c) 

3  3  3 

Zeroth  Degree  Equation 

Substitution  of  Equation  (17)  into  Equation  (15), 
omitting  again  the  higher  order  terms  of  e  yields 

N 

V2(Pq)  =  -  6(5  -  5Q)  l  c?  6(6  -  6.)  (18) 

j  =  l 

as  the  equation  of  the  zeroth  degree  in  e  where 

2  9 2  (  )  32(  ) 

V2(  )  =  - —  + - y-  .  (19) 

35^  30 

Since  this  corresponds  to  the  symmetric  case  (e  =  0, 

constant  clearance) ,  and  the  orifices  are  equally  spaced, 

the  values  of  C.  are  identical  (=  C^)  for  a  uniform  source 

3 

pressure;  and  thus  it  is  only  necessary  to  consider  the 
sector : 


2  (j  ~  1)  ir  2  j  7T 

-  i  6  i  -  • 

N  N 

Therefore  Equation  (18)  reduces  to 

V2(p2)  =  -  c°  6(5  -  50)  6(6  -  6.)  (20) 
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with  the  boundary  conditions 


po  '  1 


at 


3<po> 

H 


L 

Z  =  -  , 

D 


=  0  at  £  =  0  , 


=  P 

e.  o 

i 


where 


2i7T 

e .  =  -  ;  i=  1,2,3.  .  .,N. 

1  N 


The  solution  to  the  above  equation  was  found  in  Reference 
[15]  to  be 


2  v  l  N  L 

Pg=l  +  — <=( - €>)  + 

0  7T  1  2  D  > 


oo  cos(nN9)  cosh(nN£  )  sinh  [nN  (—  -  £  )  ]' 

l  (-Dn  - 5 - 


n=l 


n  cosh(nN  — ) 
D 


(21) 


where 


C<  =  minimum  of  (£,Cq) 
E>>  =  maximum  of  (£,5q) 


First  Degree  Equation 

The  first  degree  terms  in  e  resulting  from  combin¬ 
ing  Equations  (15)  and  (17)  yield 
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v2(p0Pi)  =  v2^  +  u2  +  u3) 


(22) 


where 


2  ^  ^2 
V  U  =  -  sin  6  —  (P~)  , 

1  2  80  U 


(22a) 


N 


V2U2  =  -  C°  cos  6  6(5  -  50)  I  6(0  -  0.)  , 
2  j=l 


(22b) 


v2u3  = 


1  n 
-  -  6(5  -  5n)  I  6(0  -  0.)  , 

9  u  J  J 

j=i 


N 


(22c) 


with  the  boundary  conditions. 


P0P1  "  ° 


at  5  =  -  , 
D 


P0P1 


0=0 


=  P0P1 


0=2  7T 


—  (PnP-,  )  =0  at  5=0. 

85 

The  solutions  of  the  above  equations  were  found  in  Reference 
[15]  to  be 


U,  = 


3NC 


4  TT 


l  (-1) 


cos  a  0 
my  m 


m=l 


|  am ( 2mN  -  1 ) 


mN  cosh  (a  5<)  sinh  am(^-5>) 


cosh (a  — ) 
m  d 


am  cosh(mN  5<)  sinh  mN  (g-  -  5>) 


cosh (mN  — ) 
D 
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cos  3  9 

_ m _ 

3  (2mN  +  1) 
m 


L 

mN  cosh  ( 3m  5<  )  sinh  ^  -  5>) 


cosh  (3_  — ) 
m  D 


3m  cosh(mN  £<)  sinh  mN  (g-  -  £>) 


cosh{mN  — ) 
D 


(23a) 


U2  =  - 


3NC 


4tt 


I  (-D 


m 


m=l 


cos(am  6)  cosh  (a^  ?<)  sinh  ^(jj-^) 


a  cosh  (a  — ) 

m  m  j-, 


l  <-l)m 

m=0 


cos(3m  9)  cosh ( 3m  5<)  sinh  3m  (jj  -  £>) 


3_  cosh(3m  -) 
m  m  j-j 


),  (23b) 


U. 


N 


=  —  l  cj 

2tt  .  .  J 
3  =  1 


1  L 

2  D 


~  (~  -  ?>)  + 


l 

m=l 


cos  (m  0.)  cos  (m  0)  cosh(m  £.)  sinh  m  (=-  -  £.  ) 

_ 1 _ ■  _ ' _ 2; _ U 2 _ 


m  cosh (m  — ) 
D 


, (23c) 


where 


a  =  mN 
m 


1  , 


3m  =  mN  +  1  . 
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If  the  source  strength,  C*3,  is  known,  Pg  can  be  cal¬ 
culated  from  Equation  (21) .  Now  if  the  perturbed  source 
strengths,  C^,  are  known,  the  perturbed  pressure,  P^,  can 
be  determined  by  summing  U^,  and  in  Equation  (23)  and 

dividing  by  Pg. 


Calculation  of  Source  Strength 

The  film  mass  flow  rate  at  the  k  source  can  be 
determined  from  the  isothermal  assumption  and  Equations  (2) 
and  (14)  to  be: 

3  3  2 

ET  —  (P  )  ds  (24) 

3n 

where  — -  is  the  outward  normal  derivative  and  the  integra- 
3n 

tion  is  performed  along  a  path  enclosing  the  kr  source. 
Substituting  the  perturbation  series.  Equation  (17),  into 
Equation  (24)  yields 


"k 


p  c 
24y  RT 


* 


mk  = 


where 


(25) 


2  3 

p  c 


24y  RT 


—  (P«)  ds  , 
k  9n  U 


(25a) 


2  3 

p  c 


24y  RT 


* 

k 


3  3 

2,  —  (P  P  )  +  3  cos  e  — 
3n  U  1  3n 


ds  . 
(25b) 
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Applying  Green's  theorem  and  Equation  (20)  to  the  integral 
in  Equation  (25a)  gives 


j>  —  (P2)  ds  =  //  V2(P2)  d6d£  =  -  C°  , 
k  9n  u  u 

and  by  Green's  theorem  and  Equation  (22),  the  integral  in 
Equation  (25b)  reduces  to 


6  2 —  (PnP,  )  +  3  cos  6  —  (P2)  ds  =  -  cf  . 

k  L  3"  01  3n  0  J  * 


Therefore  Equation  (25)  can  be  written  as 


2  3 

Pa  c  0  i 

m,  =  — - -  (CU  +  eCp  .  (26) 

24y  RT 


Expanding  the  orifice  mass  flow  in  a  Taylor's  series 


about  P 


f  “  P0 


f  and  h  =  c,  omitting  the  f  for  clarity,  gives 


3  m, 


mk(P,h)  =  mk(PQ,c)  +  (P  -  PQ) 


3P 


3m, 


+  (h  -  c). 


Vc 


3h 


P0'C 


(27) 


The  derivatives  in  the  above  equation  can  be  determined  from 
Equation  (7)  and  Equation  (8)  to  be 


3  m. 

_ k 

3P 


=  mk (po 

C 


c) 


W  , 


(27a) 
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f 0,  orifice  compensation 

Lm,  (P.;c)/c;  inherent  compensation 
K  U 


(27b) 


where 


and  the  term  in  braces  is  zero  for  choked  flow.  Since 

(P  -  P_),  is  ePi  and  (h  -  c)  is  ce(cos  8),  Equation  (27) 
uk  k 

can  be  written  as 


mk(P,h)  =  mk(PQ,c)  +  £Plk  mk(PQ/c)  >?/V+  (3-a)  e  cos  6  mk  (p  ,c)  . 

(28) 

Equating  the  constant  terms  of  Equations  (26)  and  (28)  and 
using  Equation  (7)  to  evaluate  mk(P^,c)  gives 

p0  f  N  n 

AV  G  ( - LE)  =  —  Cu  .  (29) 

V  4tt 


t  h 

The  angular  location  of  the  k  source  can  be  written  as 


it  (2k  -  1) 

8  =  -  .  (30) 

N 


Substituting  this  into  Equation  (21)  yields 
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1  + 


TT 


N  l  cosh  (nN  E,  )  sinh  nN  (—  -  E,  )~] 

-  -  g>  +  I  - 5 - 5 - 1 

L  2  D  L 

n=l  n  cosh (nN  -) 


Now  to  determine  Pq  at  the  edge  of  the  feeder  hole,  the  term 
in  brackets  is  first  evaluated  at  E,  =  £q  +  (d/D)  ;  then  it  is 
evaluated  at  ^  ^  -  (d/D);  and  the  two  values  are  averaged. 

Defining  this  value  as  Sq,  Equation  (31)  can  be  written  as 


C1^  can  be  determined  by  a  simultaneous  solution  of  Equations 
(29)  and  (32)  . 

Now  equating  the  coefficients  of  e  in  Equations  (26) 
and  (28)  and  again  evaluating  m^(PQ,c)  from  Equation  (7) 
yields 


AV  G  ( 


->  f 


¥ 


N 


*1  -  +  (3-a)  cos  6,  =  —  C, 

k  v  J  4tt 


(33) 


Introducing  Equation  (29)  into  the  above  equation  gives 


¥ 

Pi  —  +  (3-a)  cos  0. 

L  k  V  k. 


=  c: 


(34) 


,0  ... 


where  C  is  a  known  value.  Multiplying  Equation  (34)  by 
Pq ,  substituting  for  0^  from  Equation  (30),  and  rearrang¬ 
ing  yields 
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(pnpi )  = 

0  1  k 


po  v  „i 


C,t  -  (3-a) 


¥  c 


0  k 


P.  V 


cos 


t(.2k  -  1) 


N 


(35) 


(PnP.)  can  be  obtained  from  Equation  (23)  to  be 

u  1  k 


(P0PX)  =  (SJ  +  (U0 )  _  , 


3'k 


(36) 


where  (S, )  is  the  sum  of  U,  and  U„  which  have  been  evalu- 
1  k  12 

ated  at  (£q  +  d/D,  8^.)  and  (£q  -  d/D ,  0k)  ,  and  the  values 

averaged.  Substituting  7T  ( 2 j  -  1)/N  for  the  angular  coordi- 
th. 

nate  of  the  j  source,  0 ^ ,  along  with  Equation  (30)  into 
Equation  (23c)  yields 


(U_) 

3  k 


N 

I 

j=l 


T  . ,  C1  , 

3k  j  ' 


(37) 


where 


1  fl  L 

Tik  =  —  "  M 

3 K  2tt  L2  D 


oo  COS 


mir  (2j  -  1) 


N 


cos 


mir  (2k  -  1) 


N 


cosh  m  £<  sinh  m(--  £>,)_, 


m=l 


m  cosh  (m  — ) 
D 


(38) 


which  can  also  be  evaluated  at  (£q  +  d/D)  and  (£q  -  d/D)  and 
averaged.  Combining  Equations  (35)  through  (38)  yields  the 
following  system  of  N  linear  equations: 


Ajk  Cj  “  Bk  ' 


(39) 
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where 


0  if  j  ^  k 
1  if  j  =  k 


(39a) 


P  V  tt  (  2k  -  1) 

B,  =  -  (S,  )  -  (,3-a)  -  cos  -  .  (39b) 

k  1  k  4/  n 

Thus  the  values  of  can  be  determined  by  the  use  of  matrix 

J 

albegra . 

Now  with  the  source  strengths  known  the  pressures 
can  be  obtained  at  any  location.  A  numerical  area  integra¬ 
tion  can  be  performed  on  (P^  cos  0)  to  obtain  the  bearing 
stiffness,  (load/e).  Although  the  e-perturbation  analysis 
is  valid  for  small  e  only,  the  stiffness  of  journal  bear¬ 
ings  is  usually  constant  up  to  e  =  0.5. 

V.  THRUST  BEARING  ANALYSIS 

The  theoretical  analysis  of  thrust  bearings  has 
also  been  improved  during  the  past  decade  by  using  more 
exact  mathematical  models.  More  exact  solutions  of  the  gas 
film  equations  can  be  obtained  for  thrust  bearings  than 
journal  bearings  because  the  film  thickness  is  not  a  func¬ 
tion  of  the  bearing  coordinates.  An  analysis  of  a  point 
source  fed  thrust  bearing  from  Reference  [15]  will  be  used 
in  this  study.  Figure  5  shows  a  double  acting  (two  sided) 
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b.  Radial  pad. 


c.  Rectangular  pad. 
Figure  5.  (continued) 
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thrust  bearing  model  with  the  nomenclature  used  in  this 
analysis.  It  is  assumed  that  the  orifices  are  equally 
spaced;  thus,  it  is  only  necessary  to  analyze  one  radial 
segment  of  the  bearing  as  shown  in  Figure  5b. 


X  =  x/a  , 

with  boundary  conditions , 

2  b 
P  =  1  at  Z  —  ±  — 

2a 

3(P2)  1 

-  =0  at  X  =  ±  - 

3X  2 

as  shown  in  Figure  5c.  The  external  boundary  conditions  can 
be  satisfied  by  superimposing  the  Laplacian  solution  of  an 
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infinite  number  of  equidistant  (separated  by  a  distance,  a) 
sources  distributed  along  the  x  axis  and  the  Laplacian  solu¬ 
tion  of  an  infinite  number  of  equidistant  (separated  by  a 
distance,  b)  alternate  sources  and  sinks  distributed  along 
the  z  axis.  The  resulting  pressure  distribution  is 


P2  =  1 - C  7  (-l)n  Jin 

2 

n=-oo 


cosh  2tt£z  +  n  “  cos  2tt  X 

coshQr(— )  (2n  +  1)]]  -  1  J 


(41) 


where  C  is  the  source  strength. 

The  rectangular  pad  can  be  transformed  into  the 
radial  segment  shown  in  Figure  5b  using  the  transformation: 


N 

Z  =  —  Jin  R  ,  (42a) 

2  TT 


N 

X  =  —  0  ,  (42b) 

2tt 


where  N  is  the  number  of  orifices  and  R  is  the  dimensionless 
radius  (r/r  ) .  The  pressure  in  the  radial  segment  Can  thus 
be  written  as 


P2  =  1  -  -  C  l  (-l)n  Jln< 


n=-°° 


|cosh  N 

r  rll 

Jin  R  +  n  Jin  (— ) 

2  _ 

-  cos  N0 

^  cosh 

1  rl  1 

N  (n  +  -)  Jin  (— )  -  1 

2  r  2  J  ; 

(43) 


Since  the  orifice  was  located  at  the  center  of  the  rectang¬ 
ular  pad,  the  transformation  given  in  Equation  (42a) 
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evaluated  at  (Z,R)  equal  to  (b/2a,  r^/r  )  and  (-  b/2a,  r2/rc) 
results  in 


rc  = 


(44) 


An  analysis  of  a  pad  in  which  the  orifice  is  not  centered 
is  given  in  Reference  [15);  however,  it  will  not  be  con¬ 
sidered  here. 


Calculation  of  the  Source  Strength 

Since  the  clearance,  H,  is  constant,  the  mass  flow 
per  segment  can  be  determined  from  Equation  (24)  to  be 


m 


k 


24 y  RT 


6  —  (P2)  ds  . 

k  3n 


(45) 


If  the  integration  is  performed  along  the  boundary  of  the 
rectangular  pad,  the  above  equation  reduces  to 


m 


k 


12  y  RT 


r2  > 


.1 

2 


9Z 


dX  . 


Evaluating  the  derivative  in  the  preceding  equation  from 
Equation  (41)  yields 


2  ,3 

p  h  7T 

m  =  — -  C  . 

K  12y  RT 


The  clearance,  h,  can  be  written  as 


(46) 


h  =  c  ( 1  ±  e )  , 


(47) 
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where  c  is  the  clearance  in  the  centered  position,  and  e  is 
A/c  as  shown  in  Figure  5a,  page  28.  The  term,  (1  +e) ,  is 
used  for  the  larger  gap  side  and  (1  -  e)  is  used  for  the 
smaller  gap  side.  Combining  Equations  (46),  (47)  and  (7) 

(evaluated  at  the  orifice  exit)  gives 

2A  V  G  (-=£) 

- =  C  ,  (48) 

N (1  ±  e) a 

where 


a 


3;  orifice  compensation 
2;  inherent  compensation. 


Now  evaluating  Equation  (43)  at  the  edge  of  the  orifice 
gives 


2  1  C 

(P%  =1 - (S  )  ,  (49) 

2  it  u 

where  S  is  the  value  obtained  from  the  summation  term  of 
u 

Equation  (43)  evaluated  at  the  orifice  edge.  The  value  of 
C  can  be  obtained  by  a  simultaneous  solution  of  Equations 
(48)  and  (49).  A  value  of  C  is  determined  for  each  side  of 
the  bearing  and  thus  the  pressure  distribution  on  each  side 
can  be  calculated  from  Equation  (43) .  The  pressure  can  be 
numerically  integrated  to  obtain  the  load  on  each  side  and 
the  difference  gives  the  load  capacity  of  the  bearing  at 
a  given  e. 
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VI .  COMPUTER  PROGRAMS 

The  resulting  solutions  from  the  journal  and  thrust 
bearing  analyses  presented  in  the  two  preceding  sections 
involve  very  complex  equations  with  iterative  processes 
necessary  to  determine  the  source  strengths.  Thils  to  be 
of  any  practical  use,  the  equations  must  be  solved  by 
means  of  a  digital  computer.  Computer  program  listings 
obtained  from  an  IBM  704  computer  are  given  in  Reference 
[15]  for  these  analyses.  Some  modifications  were  made  to 
these  programs  and  they  were  adapted  for  use  on  the 
CDC  1604-B  computer  located  at  the  VKF .  The  modifications 
involved  correcting  several  errors  in  the  journal  bearing 
analysis  and  adding  the  capability  to  make  the  computations 
for  inherent  compensated  bearings  (i.e.,  the  last  term  in 
Equation  (27)  was  added).  The  theoretical  correction  of 
two-thirds  found  in  the  ID  analysis  to  correct  orifice  com¬ 
pensation  to  inherent  compensation  produces  the  same  load 
results  as  the  present  analysis  but  does  not  define  the 
actual  flow  obtained. 

The  program  given  in  Reference  [15]  for  the  thrust 
bearing  did  not  have  the  capability  to  compute  the  source 
strength  which  was  an  input  to  the  program.  This  was  added 
in  such  a  manner  that  both  inherent  and  orifice  compensated 
bearings  could  be  analyzed. 
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VII.  DYNAMIC  ANALYSIS 


Stability  Criteria 

Lund  [7]  shows  that  the  additional  bearing  stiffness 
obtained  from  the  "squeeze  film"  effect  is  negligible  for 
externally  pressurized  bearings  operating  at  moderate  spin 
rates  and  large  supply  pressures,  while  Pan  [17]  shows  for 
these  conditions  that  the  additional  stiffness  obtained 
from  the  "hydrodynamic"  effect  can  also  be  neglected.  Thus 
the  first  natural  frequencies  in  the  bearings  considered 
herein  can  be  calculated  using 

u>r  =  /Kr/m  ,  (50a) 

ooA  =  /KA/m  ,  (50b) 

wa  =  /Ka/I  .  (50c) 

Film  damping  data  are  also  presented  in  Reference  [7]  in 
graphical  form  such  that  one  can  estimate  if  a  bearing 
can  operate  at  (or  pass  slowly  through)  the  natural  fre¬ 
quency  with  a  given  amount  of  unbalance . 

Friction  Torque 

The  rolling  friction  (i.e.,  tare  for  roll  damping 
balance)  between  the  stationary  and  movable  parts  of  the 
bearing  can  be  determined  from  Newtons  law  of  friction, 

9u 

Tq  y  •  (5i) 

9y 
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The  velocity  distribution  in  the  8  direction  generated  from 
the  external  pressurization  is  neglected  since  this  effect 
would  be  integrated  out  during  the  following  analysis  due 
to  symmetry.  Thus  Couette  flow  is  obtained  and  the  veloc¬ 
ity  distribution  for  the  journal  and  thrust  bearings  is 

u  =  r  a)  (y/h)  ,  (52) 

where  oo  is  the  angular  velocity  of  the  rotating  part  of  the 
bearing  and  y  is  the  coordinate  normal  to  the  bearing  surface. 

Since  h  =  c  (1  +  e  cos  0)  for  the  journal  bearing, 
Equations  (51)  and  (52)  can  be  combined  to  give 

r  oo  y 

x  =  -  .  (53) 

c  (1  +  e  cos  0) 

Integrating  (r  xQ)  over  the  bearing  surface  produces  the 
angular  viscous  damping-moment  parameter, 

2it  r^  y  L 

T  =  -  (journal)  .  (54) 

2 

c  V  1  -  e 

Since  h  =  c (1  ±  e)  for  the  thrust  bearing,  combining 
Equations  (51)  and  (52)  yield 


r  oo  y  r  oo  y 

x„  =  -  +  -  .  (55) 

c  ( 1  +  e )  c(l-e) 


Integration  of  (r  Xg)  over  the  bearing  surface  gives 

4  4 

it  y(r,  -  r„) 

=  - 2 -  (thrust)  .  (56) 

c (1  -  e  ) 
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CHAPTER  III 

DESIGN  PROCEDURE 

The  gas  bearing  balance  was  to  be  contained  in  a 
space  envelope  of  2.5  inches  in  diameter  and  8.0  inches  in 
length.  It  was  required  to  support  a  radial  load  of  350 
pounds,  located  ±  2.0  inches  from  the  journal  bearing  center, 
and  an  axial  load  of  150  pounds.  It  must  also  remain  dynam¬ 
ically  stable  at  rotational  speeds  up  to  5000  revolutions 
per  minute  with  a  typical  model  installed. 

In  order  to  have  a  very  rigid  outer  race  for  the 
journal  bearing  and  to  stay  within  the  aforementioned  space 
envelope,  the  journal  bearing  inner  race  diameter,  D,  was 
selected  to  be  1.80  inches.  The  total  length  of  the  journal 
bearing  was  limited  to  6.00  inches  since  the  remaining  length 
was  needed  for  the  thrust  bearing.  With  the  space  available, 
the  maximum  thrust  bearing  size  was  an  outside  diameter,  2r-^ , 
of  2.10  inches  and  an  inside  diameter,  2r^,  1.05  inches. 
Nitrogen  was  used  as  the  lubricant. 

I.  JOURNAL  BEARING  DESIGN 

In  order  to  support  the  moment  loading,  it  was  de¬ 
cided  to  make  two  journal  bearings  of  length,  2.70  inches, 
separated  by  0.60  inch  long  undercut  to  vent  the  mass  flow 
between  the  bearings  (Figure  6) .  Since  the  journal  bearing 
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inner  race  has  a  reasonably  small  diameter  for  such  large 
loads,  the  structural  bending  had  to  be  considered.  Kenner- 
tium-W2^^"  was  the  material  selected  for  this  part  because  it 
had  a  large  modulus  of  elasticity,  54  x  10°  lb/in  was  di¬ 
mensionally  stable,  and  reasonably  easy  to  machine.  Figure 
6  shows  that  even  with  the  high  stiffness  material,  the  bend¬ 
ing  under  maximum  load  is  significant  when  considering  the 

small  clearances  required  for  gas  bearings.  The  outer  shell 

-4 

has  essentially  no  deflection,  and  therefore,  1.3  x  10 
inches  of  the  clearance  is  taken  up  by  deflection. 

Circumferential  Orifice  Spacing 

Figure  7  shows  the  effect  of  the  circumferential 
spacing  of  the  orifices  on  the  dimensionless  stiffness, 

W/e,  with  the  mass  flow  rate  held  constant.  The  load 
capacity  increases  as  the  number  of  orifices  is  increased 
up  to  about  six.  For  larger  numbers  of  orifices  the  load 
remains  constant.  Eight  orifices  per  plane  was  selected 
for  the  present  design. 

Reynolds  Number 

As  mentioned  earlier,  the  theoretical  analysis  used 
assumed  negligible  inertia  terms  and  laminar  flow.  For 
these  assumptions  to  be  true  the  Reynolds  number  must  be 
held  below  a  given  value.  In  particular,  the  inertia 

^R  indicates  registered  trademark 
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effects  can  be  ignored  if  the  modified  Reynolds  number,  Re  , 
is  much  less  than  one,  where  Re  can  be  defined  in  terms  of 
the  Reynolds  number.  Re,  to  be 

*  2h 

Re  =  Re  ( - )  ,  (57) 

L  -  L± 

with 

Re  =  p  Vm  h/y  /  (58) 

=  the  average  velocity  across  the  film. 

Gross  (18]  states  that  the  error  due  to  neglecting  inertia 
effects  is  only  about  10  percent  for  Re  =  500  and  about 
30  percent  when  Re*  =  3500. 

The  Reynolds  number  at  which  transition  from  laminar 
to  turbulent  flow  begins  is  vague.  Constantinescu  [19] 
suggests  that  Re^  =  500;  Gross  [20]  suggests  Re^.  =  1000; 
and  Shires  [21]  in  an  experimental  investigation  with  line 
source  fed,  one-dimensional  flow,  parallel  plates  found 
that  Re^  =  1300  when  h  was  .003  inch  and  Re^  increased 
slightly  as  h  decreased.  The  bearing  reported  in  Reference 
[2]  has  been  operated  with  a  double  plane  admission  inner 
race  at  a  Reynolds  number  of  1100  with  a  laminar  film.  It 
is  desired  to  stay  in  the  laminar  flow  regime  not  only  to 
be  able  to  apply  the  laminar  theory  but  to  keep  the  bearing 
friction  torque  small.  When  turbulent  flow  exists,  the 
effective  viscosity  is  increased  and  thus  the  friction 
torque  is  increased  as  shown  in  Equation  (54). 
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The  Reynolds  number  as  given  in  Equation  (58)  can 
be  expressed  in  terms  of  mass  flow  rate  as 

Re  =  M/2tt  Dy  ,  (59) 

since 


M  =  2  p  V  it  Dh  . 

Restrictor  Coefficient  Effect 

The  effect  of  the  restrictor  coefficient,  A,  on  the 
dimensionless  stiffness  is  given  in  Figure  8.  The  restric¬ 
tor  coefficient  is  essentially  a  measure  of  the  ratio  of  the 
flow  resistance  of  the  gas  film  to  the  flow  resistance  of 
the  inlet  orifices.  The  curve  shown  is  actually  for  the 
finally  selected  design.  A  family  of  similar  curves  can 
be  obtained  by  changing  the  parameters  listed  on  Figure  8. 

If  a  particular  bearing  has  a  "reasonable"  number  of  ori¬ 
fices,  the  parameters  in  Column  I  (Figure  8)  are  fixed  and 
the  abscissa  changed  to  A/V,  the  family  of  curves  practically 
collapse  into  the  curve  shown,  especially  for  values  of  V 
greater  than  five.  If  the  stiffness  is  multiplied  by  1.5, 
the  curve  also  applies  to  orifice  compensation.  Numerous 
design  charts  are  presented  in  Reference  [16]  in  this  form 
over  an  extended  range  of  A/V. 

The  characteristic  maximum  in  the  curve  is  explained 
well  in  [14] .  Briefly,  if  there  is  very  little  pressure 
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drop  through  the  orifices  (large  A)  ,  the  pressure  at  the 
inlet  to  the  gas  film  is  nearly  the  supply  pressure  and  the 
pressure  cannot  increase  significantly  when  the  bearing  is 
deflected.  When  there  is  a  large  pressure  drop  through  the 
orifices  (small  A) ,  the  pressure  level  at  the  gas  film  inlet 
is  so  small  that  very  little  load  capacity  can  be  obtained. 

For  the  present  design,  it  is  desirable  to  minimize 
the  mass  flow  rate  for  two  reasons: 

1.  Excessive  mass  flow  could  interfere  with  the 
model  motion  and  thus  introduce  unpredictable 
errors  in  the  data  obtained  from  the  balance. 

2.  In  order  to  keep  the  Reynolds  number  small,  the 
mass  flow  must  be  minimized,  Equation  (59) . 

An  examination  of  Figure  8  reveals  that  to  obtain  near 
optimum  stiffness  and  minimize  mass  flow  rate  (i.e.,  mini¬ 
mize  Nd) ,  the  film  thickness,  c,  must  be  minimized.  The 
clearance  is  limited  by  manufacturing  tolerances  and  for 
the  present  case,  by  the  surface  deflection  (Figure  6, 
page  38).  With  these  limitations,  c  was  selected  to  be 
0.001  inch. 

Determination  of  Mass  Flow  and  Supply  Pressure 

Figure  9  shows  the  effect  of  supply  pressure,  p  , 

s 

on  the  load  capacity  at  constant  mass  flow  rates  for  both 
single  plane  (L^  =  0)  and  double  plane  (L^/L  =  1/3)  admis¬ 
sion.  The  ordinate  in  this  figure  is  load  at  an  eccentric¬ 
ity  of  0.5  since  Vohr  [16]  states  that  a  journal  bearing 
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Figure  9.  Effect  of  Supply  Pressure  on  Load 
Capacity  with  Constant  Mass  Flow  Rate. 
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cannot  be  expected  to  operate  at  eccentricities  larger  than 
0.5  due  to  a  phenomena  termed  "lock-up."  It  can  be  seen  in 
Figure  9  that  the  load  increases  with  supply  pressure  up 
to  a  given  p  ,  above  which  the  load  remains  constant. 

There  is  little  effect  of  single  or  double  plane  admission 
on  maximum  load  capacity  for  the  two  cases  shown.  A  more 
detailed  discussion  of  admission  plane  location  will  be 
given  later. 

It  is  evident  from  Figure  9  that  in  order  to  get  the 
load  capacity  required,  the  mass  flow  and  thus  Reynolds 
number  must  be  large.  If  the  bearings  angular  stiffness  is 
neglected,  each  journal  bearing  would  have  to  support  some¬ 
what  over  350  pounds  since  their  centers  are  inside  the 
±  0.2  inch  specification.  It  was  the  opinion  of  the  writer, 
based  on  angular  stiffness  data  from  Reference  [16]  and 
experience  with  a  similar  configuration,  that  a  necessary 
and  sufficient  mass  flow  would  be  0.020  lb/sec.  A  supply 
pressure  of  414  psia  was  selected  since  the  load  has  prac¬ 
tically  reached  its  maximum  value,  and  larger  supply  pres¬ 
sures  would  increase  the  possibility  of  supersonic  flow 
(reduced  load)  in  the  film. 

Figure  10  shows  the  effect  of  the  orifice  plane  loca¬ 
tions  on  the  load  capacity  for  the  selected  configuration. 

It  can  readily  be  seen  that  for  constant  mass  flow,  the 
load  capacity  is  practically  unaffected  if  (L^/L)  is  less 
than  0.4  but  decreases  rapidly  as  the  orifices  are  moved 
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further  toward  the  ends.  A  value  of  one-third  was  selected 
for  L^/L  since  experience  with  a  similar  bearing  had  shown  a 
slight  increase  in  angular  stiffness  for  a  double  plane 
admission  bearing.  Reference  [16]  shows  practically  the 
same  angular  stiffness  for  the  two  configurations. 

Now  with  all  the  parameters  selected,  a  major 

problem  area  could  be  the  large  film  Reynolds  number 

(Re  =  1770) .  The  present  design  point  is  designated  "A" 

in  Figure  8,  page  43.  From  this  figure  the  writer's  first 

impression  was  that  a  turbulent  film,  if  it  occurred,  would 

essentially  increase  A  because  of  an  increase  in  y  and  a 

more  optimum  load  capacity  would  be  obtained.  The  modified 

* 

Reynolds  number,  Re  ,  from  Equation  (57)  was  found  to  be 
1.75  and  thus  the  error  due  to  neglecting  inertia  effects 
would  be  small. 

Orifice  Type 

Lund  [7]  has  shown,  theoretically,  that  pneumatic 
instability  can  be  expected  with  orifice  compensation  at 
certain  values  of  A  and  V  and  has  presented  this  informa¬ 
tion  as  stability  maps.  In  general  the  maps  show  that 
bearings  which  are  designed  for  near  optimum  stiffness  are 
subject  to  experience  pneumatic  hammer  for  values  of  V 
greater  than  about  10.  Thus  inherent  compensation  was  used 
in  the  present  design,  and  all  of  the  preceding  data  have 
been  for  this  type  of  orifice. 
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II.  THRUST  BEARING 

Restrictor  Coefficient 

Figure  11  shows  the  stiffness  versus  restrictor 
coefficient,  A,  curve  which  has  a  characteristic  optimum 
value  similar  to  the  journal  bearing.  Since  the  thrust 
bearing,  unlike  the  journal  bearing,  can  be  analyzed  for 
any  value  of  eccentricity,  the  stiffness  at  e  =  0  was 
selected  for  use  in  Figure  11. 

It  should  be  noted  that  the  curve  shown  is  for  eight 
orifices  per  side  (N  =  8)  of  the  bearing.  A  discussion  of 
the  effect  of  N  on  load  capacity  will  be  given  later. 

The  variation  of  the  load  capacity  with  eccentricity 
for  three  different  values  of  A  is  shown  in  Figure  12.  In 
all  cases  shown,  the  stiffness  (i.e.,  slope  of  the  curves) 
increases  with  e.  The  stiffness  (i.e.,  slope  of  the 
curves)  is  nearly  constant  in  all  cases  up  to  e  =  0.3,  and 
the  variation  with  eccentricity  becomes  less  as  A  gets 
nearer  the  optimum  value. 

As  was  the  case  for  the  journal  bearing,  it  is 
necessary  to  minimize  the  clearance,  c,  to  keep  the  mass 
flow  down  and  obtain  a  large  stiffness.  Thus  a  value  of 
0.0008  inches  which  is  feasible  to  fabricate  was  selected 
for  c . 
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Figure  11.  Effect  of  Restrictor  Coefficient  on 
Thrust  Bearing  Stiffness. 
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Figure  12.  Variation  of  Load  with  Eccentricity. 
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Mass  Flow 

With  the  clearance  determined,  the  load  can  be  deter¬ 
mined  as  a  function  of  mass  flow  rate  as  shown  in  Figure  13. 
The  load  is  given  for  a  value  of  eccentricity  of  0.5.  Ab¬ 
scissas  of  Reynolds  number  and  orifice  diameters  are  also 
shown.  The  Reynolds  number  can  be  expressed  as  a  function 
of  mass  flow  rate,  similar  to  Equation  (59),  to  give 

Re  =  M/8tt  rc  y  (60) 

Point  A  was  selected  as  the  design  point  which  was  con¬ 
siderably  over  the  150  pound  load  requirement.  The  writer 
had  not  had  any  experience  with  this  type  of  bearing  and 
it  was  felt  that  a  factor  of  safety  was  required  for  two 
reasons.  First,  there  was  likely  to  be  noticeable  drop  in 
supply  pressure  to  the  thrust  bearing  since  the  flow  had  to 
pass  through  some  small  channels  to  arrive  at  the  inlet 
orifice.  Second,  Vohr  [22]  showed  that  "entrance  effects" 
can  cause  considerable  loss  in  load  capacity  for  double 
acting  bearings.  These  entrance  effects  are  more  pronounced 
for  bearings  with  small  areas. 

Number  of  Orifices 

As  previously  mentioned,  the  number  of  orifices  per 
bearing  side,  N,  was  eight.  The  reason  for  this  number  was 
that  it  was  the  maximum  that  could  be  machined  without  com¬ 
plicating  the  design.  Figure  14  shows  that  the  load  capacity 
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Inherent  Compensation 
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Figure  14.  Effect  of  Number  of  Orifices  on  Thrust 
Bearing  Load  Capacity  with  Constant  Mass  Flow  Rate. 
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for  the  present  design  increases  as  N  is  increased  when 
the  mass  flow  rate  is  fixed.  There  is,  however,  only  a 
5.5  percent  increase  in  the  load  as  N  is  increased  from 
8  to  16;  and  thus  increasing  the  number  of  orifices 
above  eight  could  not  be  justified. 

III.  BEARING  DYNAMICS 

A  typical  model  that  would  be  tested  using  the  roll 

2 

damping  balance  would  have  a  mass,  m,  of  about  0.1  lb-sec  /in 

and  a  moment  of  inertia,  I,  about  its  transverse  axis  of 

2 

about  6  in-lb-sec  .  Using  these  numbers  along  with  the 
theoretical  stiffnesses  at  e  =  0  and  Equation  (50) ,  the 
first  natural  frequencies  in  the  translational  modes  can 
be  calculated  to  be 

cor  =  3290  rad/sec  =  31,400  rev/min 

=  2470  rad/sec  -  23,600  rev/min  . 

If  the  angular  stiffness  of  the  individual  bearings  is 
neglected,  the  first  natural  frequency  in  the  angular  mode 
can  be  calculated  to  be 

(jj^  =  700  rad/sec  =  6700  rev/min  . 

If  the  angular  stiffness  obtained  from  [16] , (454 , 000 
in- lb/rad  for  each  journal  and  398,000  in-lb/rad  for  the 
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thrust  bearing)  is  added  to  the  angular  stiffness  due  to  the 
separated  journal  bearings,  the  angular  natural  frequency 
is  increased  to 

to  =  772  rad/sec  =  7720  rev/min 

a 

In  any  case  the  bearing  is  well  above  the  5000  rev/min 
requirement. 
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CHAPTER  IV 

APPARATUS 

I .  GAS  BEARING 

The  gas  bearing  assembly  is  composed  of  five  basic 
components  as  shown  in  Figure  15.  Figure  15a  shows  a  photo¬ 
graph  of  the  components  while  Figure  15b  shows  the  bearing 

assembly.  As  aforementioned,  the  journal  bearing  core 

R 

material  was  Kennertium-W2  which  has  a  high  modulus  of 
elasticity.  The  outer  race,  thrust  bearing  core,  and  in- 
sert  were  fabricated  from  Ferro-Tic  ,  grade  C,  which  is 
a  corrosive  resistant  metal  with  excellent  dimensional 
stability  but  was  found  to  be  difficult  to  machine.  The 
bolt  was  made  of  stainless  steel. 

The  design  tolerances  were  not  held  very  good  dur¬ 
ing  fabrication.  Accurate  measurement  of  the  finished 
bearing  indicated  that  both  the  gap  (film  thickness)  and 
orifice  diameters  were  considerably  larger  than  the  design 
values  for  both  the  journal  and  thrust  bearings  as  shown 
in  Figure  15b.  There  was  also  a  considerable  variation  of 
the  clearances,  c,  over  the  bearing  surfaces  and  the  meas¬ 
ured  values  given  are  average.  The  finish  of  the  bearing 
surfaces  was  supposed  to  be  within  4  micro  inches  but 
actually  was  about  10  on  the  journal  bearing  surfaces. 

The  orifices  in  the  thrust  bearing  were  bellmouthed  with 
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the  larger  diameter  at  the  outer  edge  which  forms  the  re- 
stricting  area.  The  material,  Kennertium-W2  ,  used  for 
the  bearing  core  is  a  sintered  metal  and  the  metal  flaked 
off  at  the  edge  of  the  orifice  as  shown  on  the  inset  in 
Figure  15a.  An  accurate  measurement  of  the  orifice  diam¬ 
eters  was  impossible  to  obtain  and  the  value  indicated 
in  Figure  15b  was  determined  from  mass  flow  measurements 
which  will  be  discussed  later. 

II.  LOADING  EQUIPMENT 

A  photograph  of  the  bearing  with  both  axial  and 
radial  loads  being  applied  is  shown  in  Figure  16.  The 
gas  bearing  was  held  in  a  support  fixture.  Because  of  the 
large  radial  loads  required,  a  lever  arrangement  with  a 
mechanical  advantage  of  six  was  used  to  make  the  loading 
procedure  easier.  A  mass  flow  meter  incorporating  VKF 
calibrated  orifices-,  which  is  accurate  to  one  percent, 
was  installed  upstream  of  the  gas  bearing.  A  pyrometer 
was  used  in  conjunction  with  an  iron-constanton  thermo¬ 
couple  to  monitor  the  temperature  of  the  outer  race. 

R 

A  Sheffield  air  gage  was  used  to  measure  the  dis¬ 
placement  of  the  sleeve  with  respect  to  the  core  within  an 
accuracy  of  0.00002  inches  based  on  data  repeatability.  In 
order  to  get  a  good  measurement  of  the  journal  bearing  dis¬ 
placement,  a  pin  was  passed  through  one  of  the  center  vent 
holes  (see  inset  on  Figure  16)  and  attached  to  the  core. 
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The  stylus  of  the  air  gage  pickup  rode  on  the  pin  while 
the  body  of  the  pickup  was  attached  to  the  outer  race. 
Although  not  shown,  the  thrust  bearing  was  removed  and  a 
fixture  was  added  so  that  the  journal  bearing  core  was 
supported  from  both  ends  when  the  journal  displacement 
measurements  were  made.  This  was  done  in  order  to  eliminate 
bending  and  thus  obtain  more  reliable  data  that  could  be 
compared  to  theory.  Displacement  measurements  were  made 
only  for  a  center  load  applied  to  the  journal  bearing.  An 
electrical  circuit  was  used  to  determine  when  the  bearing 
surfaces  touched.  The  continuity  light  shown  in  Figure  16 
would  indicate  surface  contact. 

The  journal  bearing  could  be  loaded  at  any  axial 
location  and  the  light  observed  to  indicate  the  maximum 
load  capacity.  These  measurements  were  made  with  the 
entire  bearing  assembled  and  cantilevered  to  simulate  the 
loads  that  the  balance  would  experience  during  a  wind 
tunnel  test  program  (i.e.,  the  effects  of  bending  included). 

III.  FILM  PRESSURE  MEASURING  EQUIPMENT 

A  photograph  of  the  equipment  and  setup  used  to 
measure  the  gas  film  pressures  is  shown  in  Figure  17.  The 
same  support  structure,  mass  flow  meter,  and  temperature 
measuring  system  that  was  discussed  in  the  preceding  sec¬ 
tion  was  used.  A  sleeve  equipped  with  pressure  ports  was 
fabricated  to  fit  over  one  of  the  journal  bearings. 
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Figure  17.  Photograph  of  Gas  Film  Pressure  Measuring  System. 
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The  pressure  ports  were  spaced  every  0.375  inches  along 
the  length.  The  radial  clearance  between  the  sleeve  and 
the  core  was  0.0009  inches.  The  orifices  in  the  other 
journal  were  plugged  during  these  measurements.  A  pres¬ 
sure  gage  was  installed  in  the  place  of  the  thrust  bearing 
in  order  to  get  an  accurate  measurement  of  the  supply  pres¬ 
sure.  A  switching  valve  was  used  to  simplify  the  data 
gathering.  The  accuracy  of  the  gages  used  to  obtain  the 
pressure  measurements  is  estimated  to  within  ±  2  psi. 
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CHAPTER  V 

RESULTS  AND  DISCUSSION 

I .  LOAD  CAPACITY 

Mass  flow  rate  measurements  at  a  supply  pressure  of 
414  psia  indicated  the  thrust  bearing  mass  flow  rate  to  be 
0.032  lb/sec  which  is  the  value  predicted  from  laminar 
theory  based  on  actual  bearing  dimensions.  It  is  impor¬ 
tant  to  note  that  the  theoretical  pressure  ratio  across  the 
orifices  is  0.35  which  means  the  orifices  are  choked  and  the 
mass  flow  is  unaffected  by  the  downstream  pressure.  Thus, 
the  agreement  between  theory  and  experiment  may  be  mislead¬ 
ing.  The  measured  mass  flow  rate  is  however  considerably 
above  the  design  value  of  0.027  lb/sec  because  of  the  machin¬ 
ing  errors.  The  mass  flow  rate  per  journal  bearing  was 
measured  to  be  0.025  lb/sec  which  is  also  much  larger  than 
the  design  value  of  0.020  lb/sec.  A  theoretical  value  for 
the  actual  bearing  could  not  be  determined  since  the  diam¬ 
eters  of  the  orifices  could  not  be  accurately  determined 
as  mentioned  in  Chapter  IV,  Section  I. 

Figures  18  and  19  show  the  results  of  the  load  and 
displacement  measurements  of  the  thrust  and  journal  bear¬ 
ings  respectively.  The  journal  bearing  load  was  applied 
midway  between  the  two  bearings  with  both  ends  of  the 
core  supported.  The  load  plotted  in  Figure  19  is  one-half 


65 


Load 


AEDC-TR -69-204 


Eccentricity. 
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the  applied  load  (i.e.,  load  per  bearing).  It  can  be  seen 
from  these  figures  that  the  actual  load  for  a  given  eccen¬ 
tricity  is  about  30  percent  less  than  the  value  predicted 
from  laminar  flow  theory.  The  theory  shown  is  for  the 
actual  measured  dimensions  of  the  bearings  rather  than  the 
design  values.  The  thrust  bearing  theoretical  stiffness 
was  reduced  about  four  percent  from  the  design  value  while 
the  journal  bearing  theoretical  stiffness  was  practically 

unaffected  since  the  restrictor  coefficient  was  changed 

2 

very  little  (i.e.,  d/c  was  essentially  unchanged  from  the 
design  value) .  It  should  be  noted  that  the  thrust  bearing 
will  support  its  required  load  of  150  pounds  at  an  eccen¬ 
tricity  0.5.  Figure  18  also  shows  a  curve  for  turbulent 
theory  which  will  be  discussed  later. 

Figure  20  shows  the  results  of  loading  the  journal 
bearings  at  several  axial  stations.  The  actual  load  is 
considerably  below  the  required  value  when  applied  2.0 
inches  from  the  bearing  center.  Based  on  the  results 
shown  in  Figure  19,  the  bearings  are  operating  at  an 
eccentricity  up  to  0.7  before  contact  is  made. 

II.  PRESSURE  DISTRIBUTION 

Based  on  the  preceding  results  it  was  felt  that  the 
bearing  was  not  performing  adequately  and  an  investigation 
was  undertaken  to  determine  the  cause  of  the  reduced  load 
capacity.  Journal  bearing  gas  film  pressure  measurements 
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were  made  using  the  sleeve  described  in  Chapter  IV,  Section 
III.  First  the  pressure  just  downstream  of  the  orifices 
was  measured  and  since  the  mass  flow  rate  was  known,  the 
orifice  flow  equation.  Equation  (7),  could  be  used  to  deter¬ 
mine  the  effective  orifice  size.  The  orifice  diameter  was 
calculated  to  be  0.047  inches  as  was  previously  given. 

The  axial  pressure  distribution  for  e  =  0  and  0=0 
(i.e.,  midway  between  orifices)  is  shown  in  Figure  21  and 
reveals  the  pressures  to  be  considerably  greater  than  the 
level  predicted  by  laminar  flow  theory.  It  was  reasoned 
that  in  order  to  produce  this  increase  in  pressure,  the 
gas  film  flow  must  be  turbulent.  Figure  21  shows  that  the 
pressure  distribution  was  predicted  by  the  laminar  theory 
very  accurately  if  the  viscosity  and  thus  A  was  multiplied 
by  a  factor  of  1.56.  However,  this  increase  in  A  produces 
an  increase  in  bearing  stiffness  as  predicted  from  laminar 
theory  which  can  be  seen  from  Figure  8,  page  43.  Point  B 
represents  a  A  value  of  1.56  times  the  value  at  point  A, 
the  design  point. 

To  verify  the  existence  of  turbulent  film  flow  and 
determine  when  it  occurred,  an  experiment  was  conducted  to 
determine  the  pressure  distribution  over  a  range  of  mass 
flow  rates  and  thus  Reynolds  numbers.  The  mass  flow  rate 
was  varied  by  varying  the  orifice  diameters  and  also  vary¬ 
ing  the  supply  pressures.  The  orifices  were  filled  with 
an  epoxy  and  drilled  to  various  sizes.  Data  were  taken 
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at  supply  pressures  of  514,  414,  314,  and  214  psia.  The 
results  are  shown  in  Figure  22.  The  data  show  good  agree¬ 
ment  with  laminar  theory  up  to  a  mass  flow  rate,  m,  of 

g 

1.85  x  10  lb-sec/in  or  Reynolds  number  of  1000.  At  this 
point  the  actual  pressure  begins  to  increase  with  m  at  a 
faster  rate  than  that  predicted  by  laminar  theory.  Shires 
[21]  and  Gross  [20]  showed  that  the  gas  film  pressure  dis¬ 
tribution  for  one-dimensional,  turbulent  flow  can  be  ex¬ 
pressed  as 

~  0.133  y0,25  RT  (m')1*75  L 

(P2-1)t  =  -  3  2  -Z -  <61) 

h  pa  D 

The  above  equation  is  based  on  the  Blasius  empirical  rela¬ 
tion  for  shearing  stress  in  a  pipe, 

T  =  0.0395  Re0*25  p  V2  .  (62) 

If  Equation  (61)  is  evaluated  for  the  present  configuration 
at  the  bearing  center,  the  equation  shown  in  Figure  22  for 
turbulent  flow  is  obtained.  It  can  be  seen  that  this  equa¬ 
tion  is  an  excellent  fit  of  the  high  Reynolds  number  data 
and  thus  the  flow  must  indeed  be  turbulent  with  Re^  =  1000. 
Shires  [21]  conducted  experiments  with  various  slot  shapes 
at  several  film  thicknesses  and  found  that  varying  the  con¬ 
stant  in  Equations  (61)  and  (62)  gave  better  fits  with  the 
data  (a  variation  of  approximately  12  percent  was  required 
in  Shires  experiments) .  This  would  be  expected  since  the 
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Figure  22.  Effect  of  Film  Turbulent  Flow  on  Pq. 
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transition  Reynolds  number  is  subject  to  vary  for  different 
configurations.  It  is  therefore  somewhat  fortuitous  that 

Equatign  (61)  provides  such  an  exact  fit  of  the  data  but 

2  175 

the  fact  that  (P^  -  1)  is  proportional  to  m  is  suffi¬ 

cient  to  show  that  turbulent  flow  exists. 

As  given  in  Chapter  IV,  Section  III,  the  clearance, 
c,  obtained  with  the  pressure  sleeve  installed  is  0.0009 
inches  while  the  actual  bearing  clearance  is  0.0011  inches. 
It  is  apparent  from  Equations  (61)  and  (21)  that  the  effect 
of  changing  clearance  is  to  shift  the  curves  along  the 
ordinate.  Shires  [21]  showed  this  to  be  the  case  if  the 
clearance  was  not  changed  by  a  large  amount. 

III.  EFFECT  OF  FILM  TURBULENT  FLOW  ON  LOAD  CAPACITY 

It  would  be  very  difficult  and  is  beyond  the  scope  of 
this  study  to  develop  a  two-dimensional  turbulent  flow 
theory  to  determine  the  effects  of  turbulence  on  load 
capacity.  However,  a  one-dimensional  turbulent  flow  theory 
can  be  developed  for  the  journal  bearing  to  compare  with  l.D 
laminar  theory.  Also,  the  thrust  bearing  laminar  theory  can 
be  modified  slightly  to  give  good  approximate  results  for 
turbulent  flow. 

Journal  Bearing  ID  Turbulent  Flow  Analysis 

Equation  (61)  can  be  used  to  develop  a  ID  turbulent 
flow  theory  for  a  journal  bearing  similar  to  the  ID  laminar 
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theory  described  in  Chapter  II,  Section  I.  It  can  be  easily 
shown  that  Equation  (61)  can  be  expressed  as  a  function  of 
instead  of  m'  to  produce  the  same  result  obtained  in 
laminar  theory,  Equation  (5) .  Thus  the  load  per  unit 
length,  VI',  and  (dWVdP^)  for  turbulent  flow  can  be  ex¬ 
pressed  the  same  as  the  laminar  flow  case,  Equations  (6) 
and  (11)  respectively. 

Using  Equations  (7)  and  (61)  to  match  mass  flows 

yields 


<pf)  = 

T 


where 


1  75F  P  -c  1.75 
ApV  •  °  G(-£)  (L  -  L1)  (1 
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It  is  quite  apparent  that  Equation  (63)  is  considerably 
different  from  its  laminar  counterpart,  Equation  (9). 
The  derivative  of  P^  with  respect  to  e  can  be  obtained 
from  Equation  (63)  which  gives 


1 
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where  the  term  in  braces  is  1  for  choked  flow.  Now  Equations 
(10) ,  (11) ,  and  (64)  can  be  combined  to  produce  an  expres¬ 

sion  for  the  stiffness  of  a  turbulent  film 


dW| 

del. 


e=0 


ir  (1. 75  a-2.25) 
2  (V  -  1) 


<D1>  <D2,T>  ' 


(65) 


which  is  similar  in  form  to  Equation  (13)  .  Equation  (65) 
is  a  function  of  but  can  be  related  to  AT  or  mass  flow 
by  means  of  Equations  (61)  and  (63)  .  For  a  fixed  value  of 
P^,  the  stiffness  for  inherent  compensation  can  be  shown 
from  Equation  (65)  to  be  five-twelfths  the  stiffness  for 
orifice  compensation  (i.e.,  a  =  2  compared  to  a  =  3), 
whereas  the  ratio  was  two-thirds  for  laminar  flow. 

A  comparison  of  the  bearing  stiffness  of  the  present 
journal  configuration  as  predicted  by  laminar  and  turbu¬ 
lent  ID  theories  is  shown  as  a  function  of  mass  flow  rate 
in  Figure  23.  The  mass  flow  rate  was  changed  by  varying 
the  orifice  diameters.  It  can  be  seen  from  Figure  23  that 
the  optimum  stiffness  occurs  at  a  smaller  mass  flow  rate 
for  turbulent  flow  than  it  does  for  laminar  flow  as  would 
be  expected.  However,  the  peak  value  of  the  stiffness  for 
inherent  compensation  is  about  40  percent  smaller  for  tur¬ 
bulent  flow  than  for  laminar  flow.  It  is  interesting  to 
note  that  an  orifice  compensated  bearing  has  about  the  same 
peak  value  for  both  turbulent  and  laminar  flow  (actually 
the  peak  value  for  turbulent  flow  is  about  two  percent 
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smaller) .  The  curve  for  the  orifice  compensated  bearing 
follows  a  trend  that  the  writer  originally  reasoned  with¬ 
out  thoroughly  investigating  the  effect  of  turbulent  flow. 

The  strong  influence  of  the  mass  flow  rate  m ' ,  on  the 

z 

pressure  for  turbulent  flow  and  the  fact  that  the  orifice 
area  for  inherent  compensation  varies  with  the  clearance 
produced  a  significantly  different  trend  of  the  stiffness 
for  the  present  case. 

The  present  bearing  configuration  is  plotted  on 
Figure  23  using  an  X  but  would  be  located  at  the  point 
marked  as  +  if  the  gas  film  were  laminar.  Although 
absolute  levels  of  the  stiffness  is  not  very  meaningful 
because  ID  theory  tends  to  overestimate  stiffness,  the 
difference  in  the  levels  should  be  comparable  for  ID  and 
2D  theory.  The  stiffness  of  the  present  bearing  is  shown 
to  be  about  23  percent  less  for  turbulent  flow  as  compared 
to  laminar  flow  which  is  in  good  agreement  with  the  30  per¬ 
cent  reduction  found  by  experiment  (Figure  19,  page  67). 

Thrust  Bearing  Turbulent  Flow  Analysis 

An  approximate  analysis  for  a  turbulent  film  in 

the  thrust  bearing  can  be  obtained  by  a  simple  modifica- 

2 

tion  of  the  laminar  theory.  Since  (P  -  1)  is  propor¬ 
tional  to  mass  flow  to  the  1.75  power  for  turbulent  flow, 

the  source  strength,  C,  in  Equation  (43)  can  be  changed  to 
1  75 

KC  ’  for  turbulent  flow  where  K  is  a  proportionality 
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constant.  The  proportionality  constant,  K,  can  be  evaluated 
by  equating  the  pressure  for  laminar  and  turbulent  flow  at 
transition.  The  resulting  equation  is 


,1.75 
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os  I  (-1>"  M 


n=-co 


cosh  N 

r-, 

Jin  R  +  n  Jin  (— •-) 

L  r2  J 
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2  - 
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(66) 


where  C^_  is  the  source  strength  at  transition  and  can  be 
determined  from  the  transition  Reynolds  number  for  a 
particular  bearing  geometry.  The  source  strength,  C,  can 
be  determined  by  matching  mass  flows  as  was  done  for  the 
laminar  case.  A  fallacy  in  this  approximate  method  is  that 
the  Reynolds  number  is  not  constant  throughout  the  film 
and  it  is  necessary  to  assume  an  average  Reynolds  number 
as  described  by  Equation  (60).  However,  it  is  felt  for 
the  double  acting  bearing  that  this  is  not  a  bad  assump¬ 
tion  because  deviations  from  the  average  are  the  same  on 
each  side  of  the  bearing  and  the  effects  tend  to  cancel. 

This  method  was  used  to  predict  the  load  for  the 
present  bearing  and  is  shown  in  Figure  18,  page  66.  It 
shows  excellent  agreement  with  the  experimental  data. 

IV.  BEARING  DAMPING 

The  damping  of  the  bearing  with  no  load  was  deter¬ 
mined  experimentally  by  spinning  it  in  a  vacuum  tank  to  a 
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given  speed  and  measuring  the  decay  in  the  spin  rate  as  a 
function  of  time.  It  is  easily  shown  by  equating  moments 
to  angular  momentum  that 


T  = 
0) 


co„  -  w 

I  £n  (— - -) 

M1  ~  Ms 
t2  -  tx 


(67) 


The  damping-moment  parameter,  T  ,  was  found  to  be  invarient 

with  spin  rate  as  theoretically  predicted  in  Equations  (54) 

-5 

and  (56).  The  value  of  T^  was  determined  to  be  49.6  x  10 

in-lb-sec/rad.  The  theoretical  value  obtained  from  Equa- 

-5 

tions  (54)  and  (56)  using  laminar  viscosity  is  9.1  x  10 
in-lb-sec/rad.  If  1.56  times  the  laminar  viscosity  (the 

effective  viscosity  as  given  in  Figure  21,  page  71)  is 

-5 

used,  T^  theoretically  becomes  10.9  x  10  in-lb-sec/rad. 
Since  the  experimental  value  is  considerably  higher  than 
the  theoretical  value,  the  expulsion  of  the  fairly  high 
mass  flow  rate  must  be  producing  additional  damping. 


V.  FINAL  BEARING  DESIGN 


Design  Changes 

During  the  time  this  investigation  was  being  con¬ 
ducted,  the  gas  bearing  was  damaged  in  a  wind  tunnel  test 
and  a  new  bearing  was  to  be  fabricated.  It  was  definitely 
desirable  to  have  the  bearing  operate  in  the  laminar  flow 
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regime  since  the  effect  of  turbulence  was  to  decrease  load 
capacity  and  increase  the  bearing  damping.  In  order  for 
the  film  to  be  laminar,  the  mass  flow  rate  and  thus  orifice 
diameter,  d,  had  to  be  decreased.  Figure  23,  page  77,  shows 
that  if  only  the  orifice  diameter  of  the  journal  bearing  is 
decreased  to  the  laminar  flow  region,  there  would  actually 
be  a  slight  loss  in  stiffness  from  the  present  value.  How¬ 
ever,  if  the  clearance,  c,  is  also  decreased,  the  stiffness 
can  be  increased  as  can  be  seen  from  Figure  8,  page  43. 

This  decrease  in  clearance  was  not  desirable  originally  be¬ 
cause  of  the  structural  bending  of  the  inner  core,  but  it 
is  felt  now  that  the  bending  effect  is  not  nearly  as  detri¬ 
mental  to  load  capacity  as  the  effect  of  turbulence.  The 
net  effect  of  changing  d  and  c  was  to  give  the  final  journal 
bearing  a  slightly  larger  theoretical  stiffness  as  cal¬ 
culated  from  laminar  flow  theory. 

The  clearance  and  orifice  diameters  of  the  thrust 
bearing  were  also  decreased.  It  was  decided  to  decrease 
the  orifice  size  a  considerable  amount  to  obtain  a  large 
reduction  in  mass  flow.  The  net  effect  of  changes  in  the 
thrust  bearing  design  was  to  reduce  the  theoretical  stiff¬ 
ness  about  35  percent  since  it  was  felt  the  original  con¬ 
servative  design  was  not  necessary. 

The  only  dimensions  changed  from  the  original  bearing 

(see  Figure  15b,  page  59)  were  d  and  c.  Since  the  Ferro- 
R 

Tic  ,  grade  C,  used  for  the  original  bearing,  was  very  hard 
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to  machine  and  also  tended  to  rust  when  not  maintained  in  a 

dry  atmosphere,  the  parts  made  of  this  material  were  changed 

to  stainless  steel  for  the  new  bearing.  The  rotating  parts 

were  made  from  440  B  stainless  while  the  thrust  bearing  core 

R 

was  made  from  Armco  PH  13-8  Mo.  stainless  steel.  The 

R 

journal  inner  core  remained  Kennertium-W2  .  Accurate 
measurements  of  the  final  bearing  showed  that  all  critical 
dimensions  were  held  within  25  x  10  inches  and  the  bearing 
surface  finishes  were  within  4  micro  inches. 

Calibration  Results 

Figure  24  shows  the  load  versus  displacement  calibra¬ 
tion  for  the  final  journal  and  thrust  gas  bearings  to  be  in 
very  good  agreement  with  laminar  flow  theory.  The  mass 
flow  rates  of  the  journal  and  thtust  bearings  were  measured 
to  be  0.011  and  0.007  lb/sec,  respectively,  which  are  the 
values  predicted  by  theory.  These  mass  flow  rates  produce 
a  Reynolds  number  of  970  for  the  journal  and  390  for  the 
thrust  bearing.  It  can  be  seen  from  Figure  24  that  the 
thrust  bearing  supports  more  than  the  150  pound  design 
requirement. 

Figure  25  shows  the  results  of  loading  the  journal 
bearings  at  different  axial  locations.  The  final  bearing 
is  well  within  the  design  requirement.  A  fairing  of  the 
data  from  the  initial  bearing  is  also  shown  for  comparison. 
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Figure  25.  Variation  of  Journal  Bearing  Maximum 
Load  Capacity  with  Load  Location,  Final  Bearing. 
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The  friction  torque  of  the  final  gas  bearing  was 
measured  at  several  mass  flow  rates  by  varying  the  supply 
pressure,  pg,  to  the  bearing.  The  results  are  presented 
in  Figure  26  and  show  T^  to  increase  with  mass  flow  rate 
as  would  be  expected.  The  theoretical  value  is. in  good 
agreement  with  the  experimental  value  extrapolated  to  zero 
mass  flow  rate.  T^  for  the  initial  bearing  is  also  given 
on  Figure  26  at  the  operating  pg  of  414  psia  and  is  about 
2.5  times  greater  than  the  value  for  the  final  bearing. 

A  summary  of  the  dimensions  and  operating  characteristics 
of  the  initial  and  final  bearings  is  given  in  Table  I r 
page  87. 
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TABLE  I 


COMPARISON  OF  INITIAL  AND  FINAL  BEARINGS 


Initial 

Bearing 

Final  Bearing 

Bearing  Parameter 

Journal 

Thrust 

Journal 

Thrust 

Length,  L,  in 

2.70 

2.70 

Distance  between  orifice 

planes,  L..  ,  in 

0.90 

0.90 

Diameter,  D,  in 

1.80 

1.80 

Outer  radius,  r^ 

,  in 

1.05 

1.05 

Inner  radius,  r. 

>/  in 

0.52 

0.52 

Orifice  diameter,  d,  in 

0.047 

0.073 

0.028 

0.020 

Film  thickness. 

c ,  in 

0.0011 

0.0009 

0.00076 

0.00072 

Number  of  orifices,  N 

8 

8 

8 

8 

Temperature ,  T , 

°R 

520 

520 

520 

520 

Supply  pressure, 

■  P  ,  Psia 

414 

414 

414 

414 

Ambient  pressure,  p  ,  psia 

14.2 

14.2 

14.2 

14.2 

Restrictor  coefficient,  A 

3 .88 

9.02 

4.85 

3.86 

Mass  flow  rate. 

M,  Ib/sec 

0.025 

0.032 

0.011 

0.007 

Film  Reynolds  number.  Re 

2200 

1780 

970 

390 

Stiffness  @  e=0 , 

,  K,  lb/in 

313,000 

334,000 

755,000 

342,000 

Load  @  e=0.5,  W 

£=0.5,  lb 

170 

150 

290 

135 

Friction  torque, 

Tu,  in-lb-sec/rad 

49.6 

x  10-5 

19.9  x 

10~5 
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CHAPTER  VI 

SUMMARY 

Externally  pressurized  journal  and  thrust  gas  bear¬ 
ings  were  designed  for  use  as  a  roll  damping,  wind  tunnel 
balance.  In  order  to  be  small  in  size  and  support  the 
required  loads  necessitated  that  the  bearings  operate  at 
large  film  Reynolds  numbers.  The  first  design  resulted 
in  a  turbulent  film.  A  theoretical  and  experimental  inves¬ 
tigation  produced  the  following  results: 

1.  The  fact  that  the  square  of  the  film  pressure  is 
a  linear  function  of  the  mass  flow  rate  for  lami¬ 
nar  flow  and  is  a  function  of  mass  flow  rate  to 
the  1.75  power  for  turbulent  flow  was  verified. 

2.  Turbulent  flow  can  be  expected  at  film  Reynolds 
numbers  above  1000. 

3.  The  maximum  stiffness  for  a  given  gas  bearing 
configuration  occurs  at  a  smaller  mass  flow  rate 
for  turbulent  flow  than  for  laminar  flow. 

4.  The  maximum  stiffness  value  is  practically  the 
same  for  laminar  and  turbulent  flow  if  the  bear¬ 
ing  is  orifice  compensated  whereas  the  value  is 
considerably  reduced  for  turbulent  flow  if  the 
bearing  is  inherent  compensated. 
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5.  It  is  well  known  that  the  stiffness  ratio  of 
inherent  to  orifice  compensation  for  laminar 
flow  is  two-thirds.  It  was  shown  in  this  inves¬ 
tigation  that  the  theoretical  ratio  for  turbulent 
flow  is  five- twelfths . 

Since  it  was  necessary  to  design  the  bearings  with 
inherent  compensation  to  avoid  pneumatic  hammer  instability, 
the  resulting  turbulent  flow  of  the  original  bearing  caused 
it  to  fall  below  the  design  requirements.  A  new  bearing 
was  designed  and  fabricated  to  operate  in  the  laminar  flow 
regime.  The  results  from  this  bearing  showed  the  operating 
characteristics  to  be  much  better  than  the  original  bearing 
and  was  in  good  agreement  with  laminar  flow  theory. 
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